Coupling to optical phonons in the one-dimensional t-J model: 
Effects on superconducting fluctuations and phase separation 



S£> 

Os 



O 

m 



> 
o 

Os 

in 
o 

Os 



i 

c 

o 
o 



13 



R. Fehrenbacher* 

Materials Science Division, Argonne National Laboratory, Argonne, Illinois 604-39 

(February 1, 2008) 

The one-dimensional (ID) t-J Holstein model is studied by exact diagonalization of finite rings using 
a variational approximation for the phonon states. Due to renormalization effects induced by the 
phonons, for intermediate electron-phonon coupling, the phase separation (PS) boundary, and with 
it the region of dominating superconducting fluctuations is shifted substantially to smaller values of 
J It as compared to the pure t-J model. Superconducting correlations are weakened through charge 
density wave interactions mediated by the phonons. Possible consequences for the high-T c oxides are 
discussed. 



PACS numbers: 74.20.Mn, 74.25.Kc, 71.27.+a, 71.38.+i 



The role of the electron-phonon (e-ph) coupling for 
the mechanism leading to superconductivity (SC) in the 
high-T c cuprates is still controversial. The importance of 
electronic correlations suggests itself from the fact that 
the high-T c is achieved upon doping a Mott insulator. 
Hence, many theoretical approaches neglect the phonons 
altogether and search for a purely electronic mechanism. 
An obvious shortcoming of such theories is the fact that 
they would never lead to an isotope effect. 

While the experimental situation concerning the 
magnitude of the isotope effect and the relevant ions pro- 
ducing it is still not quite settled, it appears to be widely 
accepted by now that especially for underdoped samples, 
T c is rather susceptible especially to the mass of the oxy- 
gen ions pp. Near optimal doping, however, only a small 
isotope effect is observed. In addition, anomalous fre- 
quency shifts of phonon modes as well as other lattice 
anomalies have been observed at or slightly above T c . 
These findings strongly suggest that the electronic and 
phonon degrees of freedom should be treated on an equal 
footing. 

The simplest model to describe the electronic low- 
energy properties of doped CuC-2 sheets is the t-J model 
[3]. Due to elimination of doubly occupied sites, it 
describes strongly correlated electrons for all parame- 
ters. This explains the extreme difficulties to obtain reli- 
able results by analytical as well as numerical methods. 
An exception is ID, where an exact solution exists at 
J/t = 2 rj|. Numerically, it was possible to extract the 
full phase diagram with good accuracy 0. In brief, there 
are three distinct regions: For J/t < 2.2, the physics is 
dominated by charge density wave (CDW) fluctuations, 
followed by a small regime (2.2 < J/t < 3.2) where 
SC fluctuations are strongest, and finally, at large J/t, 
phase separation sets in. In addition, at low electron 
density, a spin gap opens in the regime of dominating 
SC fluctuations. 



In this article, we study the influence of the e-ph cou- 
pling on the phase diagram of the ID t-J model. In par- 
ticular, it is important to know how this affects the SC 
fluctuations. The e-ph coupling is incorporated by a Hol- 
stein term, i.e., & linear coupling of the local electronic 
charge to the ionic displacement. Together with the free 
phonon part, we arrive at the t-J Holstein model. Its 
phase diagram is calculated by exact diagonalization. 
The phonons are treated by a variational approximation 
which preserves their dynamics. The main results are: 
(i) Increasing the e-ph coupling shifts the regime of dom- 
inating SC fluctuations to smaller values of J/t. (ii) The 
PS boundary acquires a non-trivial density dependence 
due to an enhancement of CDW correlations by the e-ph 
coupling, and (iii) the spatial extent of SC correlations 
is shortened. 

The formation of a Zhang-Rice (ZR) singlet (repre- 
senting a Cu m oxidation state) by holes doped in the in- 
sulating half- filled state of the CuC>2 sheets is the central 
idea underlying the derivation of the t-J model from the 
more complete three-band model |3I8| . We have demon- 
strated previously in the context of CUO3 chains [H] , that 
the total energy E s of this singlet is sensitive to dis- 
placements of the planar oxygen towards the Cu atoms 
(the so-called breathing mode) through the distance de- 
pendence of the hopping matrix elements t p d and t pp . 
The dependence of E s on the displacement is partic- 
ularly strong for an in-phase motion, since the doped 
hole is predominantly in a state delocalized with equal 
probability on the four oxygens surrounding a Cu. A 
rough estimate of the coupling (as in [5]) gives a value 
of g = dE s /alQ w I0eV/A(Q the displacement). 

Experimentally, g can be estimated from the differ- 
ence 5Q of the Cu-0 bond lengths of compounds with 
pure four-fold oxygen coordinated Cu 11 (e.g., Nd2CuC>4, 
dcu-o « 1-97A E3), versus Cu m (e.g., KCu0 2 , d C u-0 « 
1.83A 11 ). From this, g = SQK « 10-15eV/A, in good 
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agreement with the theoretical estimate. K = ft 2 M as 
100eV/A 2 is the spring constant obtained from the mea- 
sured phonon frequencies [12| (M — 4Mo x , Mo x the oxy- 
gen atomic mass). 

The simplest model which incorporates the depen- 
dence E S (Q) is the t-J Holstein model [H] 
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where (i, j) denotes nearest neighbors (nn), the Cj CT , cj^ 
annihilate (create) an electron with spin er at site i, 
rii = J2 a n i<? ^ S the number, and Si the spin opera- 
tor. The projection operators 3 ,; = 1 — n^n^ elimi- 
nate configurations with doubly occupied sites, and the 
a,i,a\ are usual bosonic phonon annihilation (creation) 
operators. For simplicity, the phonon spectrum is as- 
sumed to be dispersionless, characterized by a single 
frequency ft. The e-ph coupling constant is defined as 
A = (kg 2 /2A<m) 1 / 2 . 

The total Hamiltonian has four dimensionless pa- 
rameters: the spin exchange J/t, the phonon frequency 
Ml/t, the e-ph coupling strength A/t, and the electron 
density < n < 1. The energy scale is fixed by the 
value of t. Here, we choose Ml/t = 0.2, realistic for 
the cuprates (t « 0.4eV jH], and ft » 60 - 80meV 
for the breathing modes 52)- We measure the e- 
ph coupling in terms of the polaron binding energy, 
E B = g 2 /(2K) = A 2 /(Ml) (-E B is the energy of a 
localized polaron), instead of A/t. Using the above val- 
ues for g, ft, we obtain Eb ~ 0.1 — 0.4eV, i.e., a sizeable 
E B /t w 0.25 — 1. This corresponds to a rather small 
value A = E B /W w 0.1-0.2 (W, the bandwidth in 2D). 

The model is studied using the Lanczos method on 
finite rings of length L. To minimize finite-size effects, 
we chose closed shell boundary conditions (BC), which 
prevent spurious degeneracies. For single-band ID mod- 
els, this corresponds to periodic BC if the number of 
electrons n e = 4m + 2, and antiperiodic BC if n e = 4m 
(m a positive integer). Then the ground state is a spin 
singlet. 

The infinite-dimensional phonon Hilbert space (for 
each mode) requires some further approximation even 
for arbitrarily small systems. Previous approaches to 
similar models |13II14| treated the phonons in the static 
limit, i.e., classically. The ground state within this limit 
is obtained by minimizing the electronic energy with re- 
spect to the lattice displacement (Born-Oppenheimer). 
However, neglecting the phonon dynamics can lead to a 
rather severe overestimate of the tendency towards small 
polaron formation, i.e., self-trapping of the carriers, es- 
pecially for the large optical phonon frequencies found 
in the cuprates. 



An unbiased approximation preserving the phonon 
dynamics results from a truncation of the Hilbert space 
to include n p phonon states per site. We have shown |15| . 
that for the simplest case, n p — 2, there is a natural 
and unique choice for the two states at site i, \4>i)i, |<fe)j 
defined in terms of the phonon vacuum |0) i and the co- 
herent (displaced oscillator) state |0) i = e _7?2 / 2 e' 7a * |0) i 
(n 2 = E B /{f£l)) 



\<h 
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The non-perturbative character of this ansatz allows for 
asymptotically exact ground state properties in the weak 
and strong coupling limit. For intermediate coupling, 
the approximation was checked by including additional 
phonon states (see Ref. |5J). Apart from an upward renor- 
malization of the coupling constant, we found no changes 
of ground state properties (electronic properties of the 
system with additional phonon states at E B are identical 
within an error < 0.1% to the properties of the system 
with only \4>i) l , \(j) 2 } t at E' B , with 0.55 < E' B /E B < 1 
depending on E B ). More details will be published else- 
where. Here, we shall restrict ourselves to n p — 2 using 

\4>i)i, \<h)i- 

The phase diagram of the model is determined by 
calculating (i) the compressibility k, and (ii) the crit- 
ical exponent K p which determines the decay of cor- 
relation functions in ID models belonging to either the 
Tomonaga-Luttinger liquid (TLL) |16I18| or the Luther- 
Emery liquid (LEL) universality class. Both systems 
are characterized by gapless charge excitations. The spin 
excitations are gaped in a LEL, while they remain gap- 
less in a TLL. 

The transition to a phase separated state is signaled 
by a divergent k, which we calculate from 



n 2 E {L, n e + 2) + E {L, n e - 2)- 2E {L, n e 



(3) 



Eq(L, n e ) the ground state energy of a ring with n e elec- 
trons and L sites. K p is given by K p = {TTDn 2 ^/^) 1 / 2 
(n = n e /L, the electron density) ^B], with the Drude 
weight f2Tij 



D = ttL 



d 2 E (^) 



d$2 



(4) 
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Here, Eq($) is the ground state energy for BC with a 
phase factor e 4 *, and <&o = 0, ir for periodic and antiperi- 
odic BC respectively. A value of K p > 1 indicates that 
SC correlations dominate the long-range fluctuations for 
both a TLL, and a LEL, while K p < 1 is characteristic 
for dominating CDW correlations. 
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To check the consistency of the assumed TLL scaling, 
we calculated the central charge c from the finite size cor- 
rections of the ground state energy as predicted by con- 
formal field theory, E(L)/L — £ m — irc(v c + v s )/(6L 2 ) 
|2*T] . Here v c , v s are the charge and spin velocity which 
we extract from the dynamical charge and spin structure 
factor as in Ref. |221 Using the results for L = 8, 12, we 
found a maximum deviation of 0.07 from the expected 
c = 1 for the range of parameters studied. Given the 
finite-size errors in the velocities, this confirms the va- 
lidity of TLL scaling. 

In Fig. ^ we show the calculated phase diagrams 
for E B = 0, 0.25, 0.5, 0.75, L = 10, 12, 16. Including the 
phonons, we were restricted to L < 12, but reliable re- 
sults can already be obtained for L = 10. The main ob- 
servations are: (i) With increasing e-ph coupling, the PS 
boundary, as determined from the condition 1/k = 0, is 
shifted considerably to lower J/t. An independent cri- 
terion for PS is the long wavelength divergence of the 

rijrii) 



charge structure factor S c h(k) — l/LJ2^ji e tk ^~ 1 ^' 
which gives an error estimate < 5% for the phase bound- 
ary. A similar shift without a significant reduction in 
area is observed for the region of dominating SC fluctu- 
ations occurring below the onset of PS. For Eb = 0.75, 
it extends down to values as low as J/t 1.0, compared 
to J/t w 2.2 in the pure t-J model. Qualitatively, this 
shift can be explained by an enhancement of the effec- 
tive mass of the charge carriers due to the e-ph coupling, 
i.e., a downward renormalization of t i— > t <C t, which 
results in an enhanced effective J/t. Note, that J is 
not expected to be appreciably renormalized, since the 
phonons couple to the charge. Also, the direct depen- 
dence of J on the displacements of the breathing mode 
is weak, since out of the two Cu-0 bonds between a nn 
Cu pair, they strengthen one, but weaken the other one. 
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Fig. 1. The phase diagram of the ID t-J Holstein model 
for Q,/t = 0.2 and various values of the e-ph coupling Eb- 

(ii) Especially at large Eb, the PS line and the con- 
tours of constant K p acquire a strong density depen- 
dence. In particular, a dip in the K p — 0.6 contour 
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develops for Eb — 0.75 at quarter filling. This is due to 
enhanced CDW correlations from the e-ph coupling. In 
fact, at Eb — 0.75, there is long-range CDW order for 
J/t < 0.4, as indicated by a saturation in the staggered 
density-density correlation function 
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i.e., linij^oo _D c h(j) ^ 0. Note, that for the case J/t = 
(equivalent to spinlcss fcrmions), it is fairly well estab- 
lished by analytical and numerical results |15II23| that 
the system undergoes a CDW transition at n = 0.5 
above a finite critical coupling strength Eg(Q). In 
the static limit, Eg — 0, demonstrating the impor- 
tance of quantum fluctuations. Previously |15| . we esti- 
mated E%/t w 0.53 for Q/t = 0.2, n p = 2. Note how- 
ever, that this value is definitely too low because of the 
renormalization of the coupling constant by additional 
phonon states as explained above. An analysis of the 



charge structure factor Sct(k) = l/L^j i e 



ik(j-i) 



n 3 ni) 



strongly indicates that the appearance of a 4fcp CDW is 
not unique to density n = 0.5 (which is obvious at the 
mean-field level). However, we found that for n ^ 0.5, 
the weakened lattice commensurability leads to an en- 
hanced Eg/t, as estimated from the values at which 
Sch(4fcF ) appears to diverge as a function of L. We 
conclude that the e-ph coupling in this model mediates 
an effectively repulsive interaction between the electrons, 
and hence does not support SC. 
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Fig. 2. The phase diagram of the ID t-J Holstein model 
in the Eb,J plane at quarter filling for fl/t — 0.2, and 
L = 12. 

Fig. [3 shows the phase diagram in the Eb,J plane 
at quarter filling. This plot clearly demonstrates that 
the regime of SC fluctuations is shifted almost rigidly to 
smaller values of J/t as the e-ph coupling is increased. 
Note, the small region of CDW order at small J/t, large 
Es/t. Within a strong coupling perturbation expansion 
for the case of spinless fermions, it was shown [23!, that 
apart from the well known exponential renormalization 
of the hopping, t — » te~ EB / n , the Holstein coupling also 
leads to an effective nn repulsion, which dominates in 
the limit Eb —> oo. The competition of this effective 
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phonon-mediated repulsion with the attractive spin ex- 
change term is responsible for the corresponding strong 
competition between CDW and SC correlations in the 
t-J Holstein model. For Eb > Eg (J = 0), a phase 
transition from a CDW to a metal takes place as J is 
increased from 0. 

Next, we analyze the effect of the e-ph coupling on 
the superconducting correlations directly, by looking at 
(i) the singlet SC structure factor defined as 

1 9 pair (fc) = i^e^- i ){P/P,), (6) 

3,1 

where Pj = (c/jcj+ij — c jl c j+l1^/^ creates a nn sin- 
glet pair, and (ii) the pair-pair correlation function in 
real space D pai i r (i) — (pJPq). It is important to study 
D pa i v (i) separately to distinguish short and long-range 
correlations. 

In Fig. EJa), we plot S pa i T (k = 0) at quarter filling 
as a function of J/t, Eb- We observe an overall suppres- 
sion of 5'p a i r (fc = 0), as Eb increases. This can be traced 
back to the repulsive CDW correlations introduced by 
the phonons. In addition, the peak at the onset of PS 
is increasingly smeared out as Eb grows. Looking at 
the correlations in real space, Fig. Etb), we note that 
this suppression is coming primarily from a reduction 
of long-range correlations. This underlines our conclu- 
sion that the e-ph coupling is detrimental to SC in this 
model. 




Fig. 3. (a) The singlet SC structure factor S pa i T (k — 0) 
for various Eb- The arrows mark the occurrence of PS. 
(b) The pair-pair correlation function D pa i r (i) for vari- 
ous Eb evaluated at the value J = J ma x(-Eg), where 
<S'pair(fc = 0, J, Eb) has a maximum. 

Another question concerns the existence of a spin 
gap A s , defined as the energy difference between the 
singlet ground state and the lowest triplet excitation for 
L — ► oo. It was shown [21] that the addition of a nn 
repulsion to the t-J model can lead to a large A s in the 
region of dominating SC fluctuations. Since the e-ph 
coupling also induces an effective short-range repulsion, 
one might expect an analogous appearance of a spin gap. 
Due to the restriction L < 12, an accurate finite-size 



scaling for A s is not possible. Nevertheless, extrapolat- 
ing A s (l/L) for L = 8, 12 at n = 1/4 to 1/L = 0, the 
intercept is negative or very small for all Eb- However, 
at Eb = 0.75, we consistently observe a small positive 
intercept for several J in the region of dominating SC 
fluctuations, possibly indicating a small spin gap. For 
the t-J-V model studied in we found that extrapo- 
lating results from L = 8, 12 calculations also produces a 
negative intercept as long as A s /t < 0.1. This is due to 
the fact that the scaling function is exponential at small 
1/L if a spin gap is present |2U- Assuming the same 
scaling of A s with 1/L for the t-J-V and the present 
model, we can give an upper bound A s /t < 0.1 for all 
Eb studied. 

In the 2D t-J model, coherent charge motion is al- 
ready renormalized by spin-exchange scattering result- 
ing in an effective bandwidth ~ J This should 
make the e-ph coupling more effective, i.e., probably a 
smaller coupling (as in ID) is needed to see apprecia- 
ble effects. Substantial evidence for a SC ground state 
of d x 2_ y 2 symmetry in the 2D t-J model has been re- 
ported I^SJEHI- The most reliable estimate for the SC 
phase boundary predicts a critical value J c /t sw 0.5 
above which SC appears in the interesting regime of hole 
doping (S = 1 — 71 = 0.15 — 0.25) |2H]- Our results suggest 
that the additional mass enhancement by the e-ph cou- 
pling could push down the phase boundary to lower J, so 
that SC might arise in the regime of J/t « 1/4— 1/3, rel- 
evant for the cuprates [S]. Concerning the isotope effect, 
our results are ambiguous. A smaller phonon frequency 
increases the ratio Eb/Q, i.e., strengthens the effect of 
the e-ph coupling. However, this has two competing ef- 
fects on SC: (i), the effective ratio J/t increases through 
a further mass enhancement, but (ii), the induced re- 
pulsive CDW correlations also grow. The identification 
of the dominating effect requires further calculations in 
2D. 
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